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^ . Abstract 

OO 

CO ' I study (for the first time) the dependence of Mu and the leptonic decay constant /^ on the variation 

^^ . of the running charm quark mass fhc{i^). I conclude that the present data on /^^ from Dg — > iJ,iy 

decay give a weaker constraint than Mjj, where the latter leads to the result: rhc{Mi,) — (1.08 ± 0.11) 
Qs . GeV, to two-loop accuracy in the MS scheme, in agreement with the value fhc{Mc) = (1.231q'o5) GeV 

extracted directly, within the same approximation, from Mj/^. The agreement of the corresponding 
perturbative pole mass with the one extracted directly from the data can indicate that the non- 
perturbative effects to the pole mass are negligible. Inversely, injecting the average value of rhc{Mc) 
Mh. from Af/j and Mj/^ into the ffic behaviour of fjj, I obtain /^ ~ (1.52 ±0.16)/^, which combined with 

the sum rule prediction for fajfo, gives //^^ ~ (1.75 ± 0.18)/^ in good agreement with the data. 
The extension of the analysis to the case of Jb and fo' is discussed. 
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1 Introduction 

One of the most important parameters of the standard model is the quark masses. However, contrary 
to the leptons, where the physical mass can be identified with the pole of the propagator, the quark 
masses are difHcult to define because of confinement. Some attempts have been made in order to 
define the quark pole mass within perturbation theory, where it has been shown to be IR-finite |1| 
and independent from the choice of the regularization and renormalization schemes used M . More 
recently, it has been noticed, in the limit of a large number of flavours, that the resummation of 
perturbative series can induce a non-perturbative term, which can affect the truncated perturbative 
result, and can, then, limit the accuracy of the pole mass determination [pj. However, a proper use 
of such a result needs a resummation, at the same level of accuracy (which is not often the case), 
of the Green function of a given process involving the pole mass, where some eventual cancellation 
between the resummed terms of the two series may occur. One may bypass the previous problems, 
by working, at a given order of perturbative QCD, with the running quark masses, which are treated 
like coupling constants of the QCD Lagrangian (see e.g. B), and where some non-perturbative-like 
effect is expected to be absent. A lot of efforts has been furnished in the literature 0] for extracting 
directly from the data the running masses of the light and heavy quarks using the SVZ QCD spectral 
sum rules (QSSR) Iq, 0]. In this note, I shall consider a direct extraction of the running charm quark 
mass using the observed value of Mjj — 1.865 MeV and the new data on f^)^ from Dg -^ ^v decay, 
where its average value is [g| : 

/c,~ (1.92 ±0.23)/,. (1) 

For convenience, we shall also use the fact that r = JdJId is under good control from different 
non-perturbative approaches. For definiteness and for an internal self-consistency of the analysis, we 
shall use the QSSR value §: 

^ = 1.15 ± 0.04 =^ /d ^ (1.67 ± 0.24)/,. (2) 

JD 

We shall extend the analysis to the case of fs and Jd* ■ 

2 The QCD spectral sum rules 

We shall work with the pseudoscalar two-point correlator: 

^-5(9') =1 I d^x e^«- (0|rj,(a;)jt(0)|0), (3) 



built from the heavy-light quark current: Jd{x) = {nic + 'md)c{i^Q)d, and which has the quantum 
numbers of the D meson. The corresponding Laplace transform sum rules are: 

1"°° 1 fl 

C(t)^ dt e^'^ - ImtPsit), and TZ(t) = -—log C(t), (4) 

Jt< T^ dr 

where t< is the hadronic threshold. The latter sum rule, or its slight modification, is also useful, as it 
is equal to the resonance mass squared, in the simple duality ansatz parametrization of the spectral 
function: 

- Imi/)5(i) ~ 2flM^S{t - Ml) + "QCD continuum"e(t - Q, (5) 

TT 

where the "QCD continuum comes from the discontinuity of the QCD diagrams, which is expected 
to give a good smearing of the different radial excitations Q The decay constant fjj is analogous to 
/tt — 93.3 MeV; tc is the QCD continuum threshold, which is, like the sum rule variable t, an a priori 
arbitrary parameter. In this paper, we shall impose the tc and r stability criteria for extracting our 
optimal results q The QCD expression of the correlator is known to two-loop accuracy (see e.g. 
and the explicit expressions given in pOl), in terms of the perturbative pole mass Ale, and including 
the non-perturbative condensates of dimensions less than or equal to sixn. The sum rule reads: 

l + l{^)fix)] + [cU{0,)+Ce{Oe)T]e-''>\ , (6) 

^At the optimization scale, its effect is negligible, such that a more involved parametrization is not necessary. 

^The corresponding tc value very roughly indicates the position of the next radial excitations. 

^We shall use the corrected coefficient of the quark-gluon mixed condensate given in [|ll|. This change affects only 
slightly the result. We shall also skip the negligible contribution from the dimension six four-quark and three-gluon 
condensates. Notice that there is some discrepancy on the value of the four-quark coefficient in the literature. 




with: 



X = 
C4{04) = 



9 3 

- + 2Li2(x) + logxlog(l - x) -- - log(l/x - 1) 

— log(l ~ x) + X log(l/x — 1) — {x/{l — x)) log X, 
~Mc{dd) - (a,G'2}/127r 



It can be expressed in terms of the running mass mdv) Q, through the one- loop relation 

M.(.)=^.(.){l+(f)(l + 21og^ 



Throughout this paper we shall use the values of the parameters 
used for the mixed condensate the parametrization: 

ff(rf>M-YGrd) = Mi{dd). 

We shall also use, for four active flavours [Ol ||: 

A = (325 ± 100) MeV. 



(7) 

(8) 
[2| given in Table 1. We have 

(9) 
(10) 



One can inspect that the dominant non-perturbative contribution is due to the dimension- four rric (dd) 
light quark condensate, while the other non-perturbative effects remain a small correction at the op- 
timization scale, which corresponds to r ~ 0.6 ~ 1 GcV~^ and tc — 6 ^ S GcV^. 



3 Discussions and results 
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Figure 1: Behaviour of fo versus rhc{i'). The horizontal band is the experimental domain of fo- 
The theoretical band is hmited by the two curves: (a) A=0.425 GeV, v=l GeV « r"^/^, tc > 7 
GeV^, (Jd)i/3(i GeV)=-238 MeV, M| = 0.7 GeV^, (a.G^) = 0.06 GeV^ and (b) A=0.225 GeV, 
u w Mc = 1.42 GeV, U ^ 6 GeV^ {ddy/^{l GeV)=-220 MeV, M^ = 0.6 GeV^, {a,G^) = 0.08 
GeV2. 



Given the experimental value on Md = 1.865 GeV, we present our results on fo from the first sum 
rule for different values of the charm quark running mass evaluated at p^ — v^ in Fig. 1. The second 

*It is clear that, for the non-perturbative terms which are known to leading order of perturbation theory, one can 
use either the running or the pole mass. However, we shall see that the non-perturbative effects are not important in 
the analysis, such that this distinction does not affect the result. 



sum rule gives the prediction on Mjj for each value of the charm mass (Fig. 2). The theoretical band 
is limited by the two extremal values of the QCD parameters used in Table 1. Notice that the effect 
of the errors of the different input parameters is much smaller in the ratio of sum rule 7?.(t). The 
previous analysis leads to the prediction from a two-loop calculation in the MS scheme: 



TOc(Afc) = (1.08 ±0.11) GeV, 



(11) 



where one should notice that, within the present accuracy of the data on /d^, the result comes mainly 
from the one used to reproduce the experimental value of Mu . 
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Figure 2: Behaviour of Mu versus fhdy). The horizontal band is the experimental value of M]j. The 
theoretical band is limited by the two curves: (a) the same as (a) of Fig. 1 but v—\ GeV and (b) the 
same as (b) of Fig. 1 but i/=1.42 GeV. 



We have varied the subtraction scale v in the range from r~^/^ ~ 1 GeV to M^ « 1.42 GeV. One 
should also notice, in Table 1, that the effect of t^ on the result is relatively small from the value 
ic — 6 GeV^ , where one starts to have a t stability until tc>7 GeV^ , where one has tc stability. The 
main sources of errors are due to A and v. Within the errors, the present result is in good agreement 
(though less accurate) with the value H: 



rhrXM,) = (1.23+'(J;[;t) GeV 



(12) 



obtained, within the same two-loop approximation, from Mji^. Inversely, we can use the combined 
value: 

m^lvh) ^ (1.20 ± 0.05) GeV (13) 

from Md and from Mji^ systems on the curve Jd as function of fhc given in Fig. 1. Then, one can 
deduce: 

/z,^ (1.52 ±0.16)/,, (14) 

where, as can be seen in Table 1, the errors in this determination come mainly from the perturbative 
parameters: A, rric and the subtraction scale v, and, to a lesser extent, from tc and the non-perturbative 
terms. We can compare this result with the previous value: 

/z,^ (1.35 ±0.07)/,, (15) 

obtained by using the perturbative pole mass of the charm quark propagator ||9|, |lO|. The good 



Table 1: Different sources of errors in the estimate of fo 
Sources |A(/d//^) 



A = (325 ± 100) MeV 0.12 

z^ = (1.20 ±0.22) GeV 0.08 

fhciMc) = (1.20 ± 0.05) GeV 0.05 

tc = (6.5 ± 0.5) GeV2 0.04 

(Jd)i/3(i GeV)=-(229 ± 9) MeV 0.02 

Mo^ = (0.8±0.1) GeV2 0.02 

(a^G^) = (0.07 ± 0.01) GeV^ 0.01 

Total 0.16 



agreement between the two results in Eqs (|lj) and (|lj) within the errors may be an indirect indication 
that the pole mass defined at a given order of perturbation theory, can provide a good description of 
the physical process in this channel, and that the eventual non-perturbative power corrections induced 
by the resummation of the perturbative series remain small corrections. This observation can also be 
supported by the agreement of the value of the perturbative pole mass obtained here and of the one 
from the J/iJj systems, where both values have been obtained at two-loop accuracy. A further support 
of this argument can also be provided by the agreement of the pole mass deduced from Eq. (8) using 
the value of the running mass, with the one extracted directly in p4| . Finally, using the previous 
value of the ratio fos/fo given in Eq. (2), one obtains: 

/i3,~ (1.75 ±0.18)/., (16) 

which is in good agreement with the data quoted in Eq. (1). A natural improvement of the analysis 
done in this paper is a much more precise measurement of fo^ or/and an evaluation of the QCD 
two-point correlator to three-loop accuracy. The two projects are expected to be feasible in the near 
future. 

4 Extension of the analysis to /b and f^* 

The extension of the previous analysis to the case of the B meson is not very conclusive due to the 
strong dependence of the result on the value of the subtraction scale v, which can vary, for the B 
meson case, in a larger range. Instead, using our previous observation on the agreement of the results 
on fu from the running and perturbative pole masses, we also expect that the result obtained from 
the pole mass of the b quark ^ Q should also be reliable, which is: 

/s^ (1.49 ±0.08)/.. (17) 

We have also extended the previous analysis to the case of the vector current having the quantum 
number of the D* . In this case, and working with the correlator having the same dimension as ^'5(9^) 
n, we do not obtain a r stability as the coefficients of the chiral condensates are opposite of the 
pseudoscalar ones. Then, we conclude that, from this quantity, we cannot have a good determination 

oifD'. 
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